YITP-SB-12-40 
ICCUB-12-356 



Dark Radiation Confronting LHC in Z' 

Models 



A. Solaguren-Beascoa b , M. C. Gonzalez-Garcia 



a,b,c 



& lnstitucio Catalana de Recerca i Estudis Avangats (ICREA) 

b Departament d'Estructura i Constituents de la Materia and ICC-UB 
Universitat de Barcelona, Diagonal 647, E- 08028 Barcelona, Spain. 

C C.N. Yang Institute for Theoretical Physics 
State University of New York at Stony Brook 
Stony Brook, NY 11794-3840, USA. 



Abstract 

Recent cosmological data favour additional relativistic degrees of freedom beyond 
the three active neutrinos and photons, often referred to as "dark radiation". Exten- 
sions of the SM involving TeV-scale Z' gauge bosons generically contain superweakly 
interacting light right-handed neutrinos which can constitute this dark radiation. 
In this letter we confront the requirement on the parameters of the E§ Z' mod- 
els to account for the present evidence of dark radiation with the already existing 
constraints from searches for new neutral gauge bosons at LHC 7. 



Additional heavy neutral Z' gauge bosons are predicted in many extensions 
of the standard model (SM) such as most Grand Unified Theories (GUT) and 
superstrings [lj. Generically a model with an extra U(l)' gauge symmetry 
is characterized by the mass of the Z', the £7(1)' gauge coupling and the 
chiral charges of the matter fields (in some cases including additional particles 
with exotic SM charges to cancel anomalies), and the possible mixing angle 
between the Z and Z'. Particularly well motivated and consistently anomaly- 
free constructions are Z' models based on the E 6 GUT group and for this 
reason they have been extensively studied in the literature. In what follows 
we will concentrate on these type of models though the results here presented 
can be easily generalized to other assignments of Z' charges and couplings. 

Since E 6 is a rank 6 group, it contains in general two neutral gauge bosons 
beyond those of the SM. These couple to two new hypercharges if) and \ 
corresponding to the U(l) symmetries in E 6 /SO(10) and SO(10)/SU(5) re- 
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Table 1 

The SM quantum numbers and hypercharges of the U(l) x and the U(l)^ of the 
relevant fields. 

spectively. These hypercharge quantum numbers for the SM fermions and 
right-handed neutrinos are given in Table [TJ Here we focus on the case where 
the gauge sector contains only one additional U(l) symmetry at low energies. 
Therefore there is a continuum of possible models where the new gauge boson 
couples to one linear combination of Y x and parametrized by a mixing 
angle f3 0. 

Yp = cos /3Y X + sin /3 Y^ , (1) 

In what follows will make our study for an arbitrary value of (3 which cor- 
responds to the most general single Z' model where the new £7(1) can be 
embedded in a primordial E$ symmetry. We can chose < (5 < tt since the 
charges merely change sign for /3 — > (5 + 7T. In this case the sign of the mix- 
ing angle between Z and Z' becomes physical. However the latest analysis 
of precision electroweak data jl] constraints the Z-Z' mixing angle to be at 
most 0(1O -3 ) for Mz> < 1 TeV and theoretically the mixing angle is expected 
to decrease inversely proportional to M§,. For the sake of simplicity we will 
neglect the small effects associated with the Z Z' mixing and set it to zero 
in our calculations. Finally the value of the U(l)/3 coupling is fixed by the 

condition of coupling constant unification gy^ = \J\g2 where gi is the SU(2) 
coupling constant. 

If present at the TeV scale these Z' bosons would lead to unmistakable sig- 
natures at colliders with the strongest constraints arising from searches of 
resonances decaying into dilepton final states at LHC7 [2113] • In Figured] we 
show the strongest present bounds on <j(pp — > Z') x Br(Z' — > l + l) at 95% 
CL from these searches obtained with 5fb _1 luminosity [2] together with the 
theoretical expectations in the E 6 Z' models for different values of 0. The 
theoretical expectations were obtained introducing the E 6 models in the pack- 

1 A special case that is often considered in the literature is U(l) rn which in our 
convention corresponds to /3 = tan _1 (— \/5/3) ~ 0.717T 



2 




Fig. 1. 95% upper limits on the productino ratio cr(pp — > Z') x Br(Z' — > from the 
combined di-muon and di-electron searches in Ref. [2] as a function of the resonance 
mass. Also shown in the figure the predictions for ~Eq Z' models for different values 
of the model parameter (3 as labeled in the figure. The curves corresponds from left 
to right to the values of f3 listed from top to bottom. 

age MADEVENT [5] and with the CTEQ6L parton distribution functions [6] . 
They are in perfect agreement with the corresponding theoretical curves shown 
by the collaborations for some specific values of /3. From the figure we read 
the most up-to-date bound on Mz> as a function of the model parameter /3 
which we will use in the following and it is shown in FigJH 

A TeV-scale Z' also has important cosmological implications. Since the right- 
handed neutrinos carry a non-zero U(l)' charge, the U(l)' gauge symmetry 
prevents the large right-handed Majorana masses needed for the ordinary neu- 
trino seesaw mechanism. In this case right-handed neutrinos remain massless 
while the left-handed SM states must be light to account for the oscillation 
data [7]. Alternatively right-handed and left-handed neutrinos can combine 
to form three light Dirac neutrinos. Either way the lightness of the observed 
neutrino masses is generated by some additional mechanism and the three 
right-handed states are either light or massless and only couple to the SM sec- 
tor via the Z' interactions. In this case the theory contains the right-handed 
neutrinos relativistic degrees of freedom in addition to photons and the three 
left-handed neutrinos of the SM which can contribute to the expansion rate 
of the Universe as a new form of "dark radiation" and consequently affect the 
cosmological observations. 

In particular faster expansion would lead to an earlier freeze-out of the neutron 
to proton ratio and would lead to a higher 4 He abundance generated during 



3 



Big Bang Nucleosynthesis (BBN) [8]. At later times dark radiation would also 
alter the time for the matter-radiation equality with the corresponding impact 
in the observed cosmic microwave background (CMB) anysotropies as well as 
affect the large scale structure (LSS) distributions . 

In Refs. [9j[T0] it was first discussed in the context of BBN the implications 
of a superweakly interacting light particle, such as the right-handed neutrinos 
coupling to a heavy Z' . Because of their superweak interactions, such particles 
decouple earlier than ordinary neutrinos and consequently their contribution 
to the energy density budget of the Universe - and therefore to its expansion 
- is suppressed with respect to that of the left-handed neutrinos. The cosmic 
radiation content is usually expressed in terms of the effective number of ther- 
mally excited neutrino species, iV e fj with its standard value being iV e fj = 3.046. 
After their decoupling (for T < T^ c < T^ c ) the three superweakly interacting 
light right-handed neutrinos contribute to AN c s as: 

— (t;) 4 = 3 < 2 > 

where g(T) is the effective number of degrees of freedom at temperature T. 
Neglecting finite mass corrections, g(T) = gsiT) + £g F (T), where g BF (T) 
are the number of bosonic and fermionic relativistic degrees of freedom in 
equilibrium at temperature T. Thus at T^ c ~ 3 MeV, g(T^ c ) = 43/4 cor- 
responding to three active neutrinos, e ± and photons. In calculating g(T) at 
higher temperatures one must take into account the QCD phase transition at 
temperature T c . Above T c the quarks and the gluons are the relevant hadronic 
degrees of freedom, while below T c they are replaced by the hadronic degrees 
of freedom. At present the estimated value of T c = 154 ± 9 MeV, [llj and 
the evolution of the energy and entropy across the QCD phase transition are 
obtained by means of state of the art lattice QCD simulations [12]. Using the 
results in Ref. [12] and including finite mass effects we obtain the T depen- 
dence of g(T) (without the right-handed neutrino contribution) shown in the 
left panel in Fig. [2J 

In the Z' models here discussed right-handed neutrinos are kept in equilibrium 
by their interactions with the SM fermions mediated by the Z 1 so 
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Fig. 2. Left: The effective number of degrees of freedom as a function of the tem- 
perature without including the contribution of the right-handed neutrinos. Right: 
Decoupling temperature of the right-handed neutrinos as a function of the Z' mass 
for different values of the model parameter /3 as labeled in the figure. The curves 
corresponds from left to right to the values of j3 listed from top to bottom. 

gv R =2 and n UR is the number density of a single generation of right-handed 
neutrinos which follows a Fermi-Dirac distribution f VR (k) = ek ^ +1 - s = (p + 
q) 2 = 2pk(l — cos8) is the COM energy in the collision and v — (1 — cos#) 
where 9 is the relative angle of the colliding right-handed neutrinos. 

In Eq. ([3]) the annihilation cross section is 



cr f (s) = a(u R u R //) 



n! 



c 16tt Mf , 



(^) 2 {(l + f) 



(4) 



r /3 



where is the number of colours of the fermion / and (3j = y 1 

The UrS will decouple when their interaction rate drops below the expansion 
rate of the Universe H(T). Right before decoupling 



H(T) = \ 



A^G N (g(T) + S) 



45 



4 J rp2 



(5) 



— contribution to the expansion due to the 3 
The decoupling temperature is denned by the 
condition T UR (T^ C ) = H(T^ C ). Its value as a function of M Z i and can be 
found by numerically solving this equality with the expressions in Eqs. (J3]) 
and ([5]) and g{T) as given in FigfSl We show it in the right panel of Figj2j The 



where we have included the 
generations of massless u' R s 



coupling Yp R vanishes for tan/3 = vl5, ie (3 = 0.41967T. Thus for this value 
of the z/r's are never in equilibrium with the SM particles. Consequently, 
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Fig. 3. The equivalent number of extra neutrinos AiV e g due to 3 generation of 
massless vrs as a function of the Z' mass for several values of the model parameter 
/3 as labeled in the figure. The curves corresponds from left to right to the values of 
(3 listed from top to bottom. 

as seen in the figure, — > oo independently of Mz> as /3 — > tan _1 (Vl~5). 
Conversely Yg R is maximum for = 0, ir ie for the U(l) x hypercharge so for 
a given M z > T^ c the lowest for this model. 

Once we know as a function of the model parameters we can derive the 
contribution of the right-handed neutrinos to the dark radiation at tempera- 
tures below their decoupling as Eq. (J2]). This is illustrated in Fig. [3] where we 
plot AiV e fj as a function of the Z' mass for several values of the model param- 
eter (3. As expected as (3 approaches the decoupled model (3 = tan" 1 (a/15) ~ 
0.42, the additional number of effective neutrinos becomes small (asympoti- 
cally zero) for any value of Mz'- We note that for the range of Mz> masses 
plotted in the figure, the minimum value of AN e s ~ 0.2 corresponds to the 
plateau around g(T^ c ) ~ 86 for the corresponding decoupling temperatures 3 
GeV< Tfoc < 10 GeV as can be seen in Fig. [2 

All this implies that the cosmological information on AiV e fT can be used to 
constraint the Z' properties. In particular during most of the last two decades 
the measured primordial abundances of 4 He and other nuclei were used with 
this purpose [131415] . Since at the time the observed abundances were compat- 
ible with the presence of no dark radiation, generically lower bounds on the 
Z' mass were derived. 

This situation is altered at present as recent analysis of cosmological data 
suggest a trend towards the existence of "dark radiation" (see for exam- 
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pie [T61422] ). The Wilkinson Microwave Anisotropy Probe (WMAP) collab- 
oration found N eS = 4.34±g;|| based on their 7-year data release and addi- 
tional Large Scale Structure (LSS) data [H] at la in a ACDM cosmology. In 
Ref. [IT] N eS = 4.35±J;54 was found in a global analysis including the data 
from cosmic microwave background (CMB) experiments (in particular the 
from WMAP-7), the Hubble constant HO measurement [23], the high-redshift 
Type-I supernovae [24] and the LSS results from the Sloan Digital Sky Survey 
(SDSS) data release 7 (DR7) halo power spectrum [25] in generalized cos- 
mologies which depart from ACDM models by allowing not only the presence 
of dark radiation but also dark energy with equation of state with u ^ — 1, 
neutrino masses, and non-vanishing curvature. More recent measurements of 
the CMB anisotropy on smaller scales by the Atacama Cosmology Telescope 
(ACT) [26] and South Pole Telescope [27] experiments seem to also favour a 
value of N e g higher than predicted in SM. Cosmological constraints from BBN 
indicate, as well, that the relatively high 4 He abundance can be interpreted 
in terms of additional radiation during the BBN epoch [28H30] (see [31] for a 
recent review). In Ref. [31] using observed D and 4 He abundances and fitting 
simultaneously N eS and the the baryon asymmetry, it is found N eS = 3.71^q;45 
at lo". 

This positive evidence of a non-zero AN e g can be interpreted in the Z' models 
as evidence of a TeV scale Z' and used to quantitatively derive the required 
values of Mz> as a function of the model parameter /3 and to compare those 
with the present LHC7 constraints 0- We find that the 95%CL LHC7 lower 
bounds on Mz> as derived from Fig. [1] imply that within these scenarios the 
effective number of neutrinos is bounded to be GO 
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We also show the results of this exercise in FigJD where we plot the values of 
Mz> required to generate two illustrative la ranges of N eS , N e g = 4.35^o;54 
(0.76 < AiVeff < 2.7) as obtained in the analysis of Ref. [17] (right panel) 
and for the lower range presently favoured by BBN nucleosynthesis N e g = 
3-71±a45 (0-21 < AiV cff < 1.13) (left panel). Also shown in the figure are 95% 
CL the present bound from LHC7 as obtained from Fig. [TJ As seen in the 
figure already with the existing bounds from LHC7 most of the values of M' z 
required to account for the larger amount of dark radiation (N e g = 4. 35jzJ;f 4 ) 
are disfavoured. 



2 Along these lines Ref. [32] presents an estimate of the corresponding gauge boson 
mass range for an specific model with two additional Z' bosons. See also Ref. |33j . 

3 In particular for the 77 model AN^ ax = 0.42. 
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Fig. 4. Required values of Mz< for two illustrative la ranges of N e g, N c g = 4.35_ ' 54 
(0.76 < AiV e ff < 2.7) as obtained in the analysis of Ref. p2] (left panel) and 
for the lower range presently favoured by BBN nucleosynthesis N e g = 3.71^'^ 
(0.21 < AiVeff < 1.13) (right panel). The solid red curve are the 95% CL lower 
bounds on Mz> from searches for Z' resonance in di-lepton production at LHC4 in 
Ref. [2]. 



In summary in this letter we have explored the possibility of accounting for 
the increasing evidence of dark radiation found in recent analysis of CMB, 
LSS, and BBN data in the framework of E6 Z' models in the light of the 
present LHC7 constraints on these scenarios (which we display in Fig{T]). In 
these models additional radiation exists in the form of three generations of 
relativistic right-handed neutrinos in an amount that depends on how long 
they are kept in equilibrium with the SM particles by the Z' interactions as 
we quantify in Figs. [2] and [31 We conclude that within the present 95% CL 
bounds already imposed by LHC7, these scenarios cannot account for extra 
radiation in excess of 1.25 effective neutrino species for any value of and 
in excess of 0.5 effective neutrinos for models with 0.1 < f3/n < 0.75 (see 
Eq. ([2l)). The potential of the future LHC14 runs to further constraints the 
amount of dark radiation in these scenarios can be easily read from Figj3j For 
example, should no resonance be found at LHC14, it is foreseeable that the 
lower bound 95% of Mz> will become at least ~ 5 TeV |34j . In this case the 
amount of dark radiation in these Z' models will be limited to at most that 
of 0.3 effective neutrinos. 
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